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Abstract. We consider a disordered quantum metamaterial formed by an array of 
superconducting flux qubits coupled to microwave photons in a cavity. We map the 
system on the Tavis-Cummings model accounting for the disorder in frequencies of 
the qubits. The complex transmittance is calculated with the parameters taken from 
state-of-the-art experiments. We demonstrate that photon phase shift measurements 
allow to distinguish individual resonances in the metamaterial with up to 100 qubits, 
in spite of the decoherence spectral width being remarkably larger than the effective 
coupling constant. Our simulations are in agreement with the results of the recently 
reported experiment [T]. 
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1. Introduction 

The novel type of quantum metamaterials ElE] composed of arrays of superconducting 
two-level systems (qubits) offers a platform for quantum simulators and quantum 
memories HE]. Quantum metamaterials can be employed as test bench for studies 
of fundamental phenomena as ensemble quantum electrodynamics and spin resonance 
physics on macroscopic level BIZI- The flux qubit iE], which can be viewed as a 
artihcial atom [10], is a /im-sized superconducting loop with several Josephson junctions 
acting as nonlinear circuit elements. The strongly anharmonic potential of the flux 
qubit results in an effective two-level structure of the lowest pair of energy levels of the 
system, where the excitation frequencies fall into GHz range. Ground and excited levels 
of the qubit correspond to quantum superpositions of states with opposite directions of 
macroscopic persistent currents in the qubit loop. When coupled to a photon held in a 
superconducting cavity, the qubit becomes ’’dressed” with photons. The experimental 
studies of superconducting qubit-cavity quantum electrodynamics have shown vacuum 
Rabi oscillations mm, spin echo and Ramsey fringes m m, emission of single 
microwave photons [T3], and Lamb shift |16j . 

While the artihcially made ensembles of superconducting qubits have an 
unavoidable spread of excitation frequencies, in contrast to identical natural atoms, 
they are easily tunable by varying an external magnetic hux. This fact makes possible 
observation of the fundamental phenomena as dynamical Gasimir effect by applying of 
a non-stationary external drive [13 m, under which GHz photons are created from 
cavity vacuum huctuations. In case of the multi-qubit system, a cm-long cavity is 
used as a transmission line where the virtual photon exchange provides a long-range 
qubit-qubit interaction in a sub-wavelength metamaterial. It was shown by Tavis and 
Gummings [19] that N identical two-level systems coupled to the single photon mode 
generate a collective enhancement of the coupling constant proportional to \/N. In a 
qubit array with disordered values of excitation frequencies photons are also coupled 
to a collective superradiant mode and decoupled from other N — 1 subradiant modes. 
Gollective enhancement combined with a possibility of tuning the parameters of the 
metamaterial by external magnetic field opens new possibilities for quantum information 
technology. For example, it has been suggested, that by applying gradients of magnetic 
held along the array of qubits, the system can be operated as a quantum memory with 
information encoded in collective qubit states [2U] . 

The present work is motivated by recent experiments by Macha et ah [1] 
where the collective states of an array of superconducting hux qubits were probed 
by microwave photons transmitted through the resonator, in the regime of relatively 
strong decoherence and disorder in qubit excitation frequencies. The experiment has 
shown that resonant photon phase shift ehect is quite prominent and stable against the 
decoherence. We perform simulations of the phase shift as a function of external hux 
that modihes the qubit excitation energies, in the setting similar to Ref. [T]. In our 
approach, we perform a diagonalization of the Hamiltonian in a single excitation basis. 
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We analyze the phase shift portrait for different nnmber of qnbits and varions valnes of 
disorder in qnbit excitation freqnencies. We show that photon phase shift measnrements 
are capable to resolve individnal resonances for np to = 100 qnbits. Onr simnlations 
closely reprodnce the experimental observations. 

2. Model 

2.1. Hamiltonian of Disordered Metamaterial 

Considering the interaction of N qnbits with a photon in a resonator (see Fig. [^, 
we start from the Tavis-Cnmmings Hamiltonian, which inclndes qnbit-qnbit interaction 
term 

N N N-l 

H = ujra^a + eiof a~ + Y 9i{^t « + ) + Y 9qq,i^t+iH (1) 

1=1 2=1 2=1 

Bosonic operators a'^,a in (|^ correspond to the photon mode Ur in the cavity and 
the single i-th qnbit is described in terms of nppering and lowering Pauli operators 
— k)(5'IWi~ = |5')(^I acting on its ground and excited states. The first and the 
second terms in Eq. ([^ correspond to photon and qubit excitations of frequencies Ur 
and Cj, the third one is qubit-cavity coupling written in rotating wave approximation. 
The last term is an effective direct nearest neighbor qubit-qubit interaction. 



OQDOOOOQO fluxqubits 


microwave resonator 


ooooooooo 



Figure 1: Sketch of quantum metamaterial formed by an array of flux qubits embedded 
into a microwave resonator. 

The qubit excitation frequency e follows from microscopic Josephson and charging 
energies, Ej and Ec, which imply Ej > Eq in the operating regime of flux qubit. When 
the external magnetic flux threading the qubit loop is equal to half of the magnetic flux 
quantum $ = $ 0 / 2 , with $0 = h/(2e), the total energy of the qubit has two symmetric 
minima related to the opposite circulations of persistent currents |1) and H) where the 
tunneling rate A between two wells depends on Ej and Eq. The condition that A is 






4 


higher than dephasing rate allows for qnantnm snperpositions of |J,) and |t) resnlting 
in two non-degenerate states \g), \e) being the basis of Q. By detnning the magnetic 
flux from $ = <ho/2 one shifts the two minima by the energy 


£i(®) 



( 2 ) 


where Ip^i is the z-th qubit nominal current [8]. The excitation frequency of z-th qubit 
is given by the Floquet relation 

£,(4) = v^A? + £,(4)2. (3) 

The qubit-cavity coupling constants in Eq. ([^ are renormalized bare constants 
Qi = written in the rotated basis \g), \e). Everywhere below we work with 

the effective constants. Under uniform flux biasing conditions, the spread in depends 
mainly on qubit excitation gaps Aj, rather than Jp j, due to exponential dependence of 
Aj oc y/Ej /Ec exp(—a Ej / Ec) on the Ej/Ec ratio which fluctuates from one flux 
qubit to another. The system under consideration [T] is strongly disordered because 
the spread of excitation gaps Aj is estimated as high as cta < 20%. The system has 
an intermediate collective coupling strength value hi, which is smaller than decoherence 
rate T^ but larger than relaxation rate k of the cavity 


T^ > f] > K, ^ = J 

y 1=1...N 

In the experiment by Macha et ah |1], the external magnetic flux $ tunes qubit 
excitation frequencies in resonance with the cavity mode providing the shift of the phase 
<p(<h) of a weak external probe signal at cavity mode frequency Urj^n = 7.78 GHz. From 
the measured <p(‘F) it was found that number of qubits in an ensemble, which form a 
collective mode, corresponds to A = 8. Relevant parameters of the studied metamaterial 
are the following: qubit-cavity coupling strength gi/2Ti = (1.2 ± 0.1) MHz, dephasing 
F(p/27r = 55 MHz, persistent current Ip = (74±1) nA and average value of the excitation 
gap A/27r = 5.6 GHz. The qubit-qubit interaction, estimated as gqq^i/27i < 1 MHz, has 
negligible effect on </3(‘F). The interaction only slightly shifts energies of collective states 
and does not contribute directly to the collective qubit-cavity coupling strength G, which 
is the most relevant for the phase shift. 


2.2. Exact Diagonalization Approach 

The qubit dephasing in the system under consideration leads to smearing of the photon 
density of states which becomes not informative. However, the experimentally measured 
phase locking effect is quite prominent [1]. In our work here, we focus on the influence 
of relatively large diagonal disorder in the qubit excitation frequencies and their number 
N on the phase shift (^(4*) of the transmitted photons. We calculate the phase shift 
from a complex phase of the photon Green function 93 ( 4 *) = argP^j, where u is the 
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probe frequency and $ is the external flux, counted from the symmetry point value 
$ 0 / 2 , which shifts all excitation energies in the qubit ensemble. 

Considering approximately resonant regime where the excitation frequencies of 
qubits and resonator mode are close to each other 

|ej — UJr\ < (Jr (4) 

we hnd a spectrum of the system in a regime of single excitation by means of 
exact diagonalization of the Hamiltonian. This regime is realized under experimental 
conditions due to small microwave probe power and low temperature of the system. 
Our solution corresponds to the fully quantum regime where all the variables a~, a are 
quantum helds. 

Below we introduce the basis of + 1 states of approximately equal energies which 
are related to single excitation either in the photon cavity or in one of N qubits. Namely, 
the state where a single photon is excited and all qubits are in the ground state reads 

| 1 ) = | 1 ; 0 ,..., 0 ) 

Af+l 

The excitation in 1-st, Tth or iV-th qubits without the photon are given by 

| 2 ) = | 0 ; 1 , 0 ,..., 0 ), 


|i + l) = |0; 0 ,..., 1 ,..., 0 ) 

and 

|N + 1) = |0; 0,...,0,1). 


The Hamiltonian (j^ in a matrix form 
basis |i) reads 


(Jr 

9i 

92 

9i 

ei 

9qq,l 

92 

9qq,l 

^2 

93 

0 

9qq,2 

9a 

0 

0 

95 

0 

0 


= (i|if|j) in terms of this single excitation 

93 Qa 95 

0 0 0 ... 

9qq,2 0 0 

^3 9qq,3 0 .... 

9qq,3 ^4 Qqq^^ 

0 9qq,A ^5 
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The photon Green function is given by the first diagonal {1,1} element of the 
resolvent matrix G 

Gjj = - diag(m, iT^, (5) 

This is the matrix of retarded Green functions, where the last term diag(m, iT^,...) 
introduced in order to take into account the damping in the cavity k/2tt = 715 kHz and 
dephasing of qubits. We note that the above expression for the photon Green function 
and Gi^j allowing us to calculate the phase shift coincides with the generic equations 
obtained earlier by Volkov and Fistul [3]. 

3. Results 

3.1. Strong Disorder and Large N 

In our numerical studies, we model the system with large amount of qubits and in 
the presence of disorder between them. Namely, we fix metamaterial parameters and 
plot the phase shift increasing qubit number N, but keeping spread of Aj constant, 
see Fig. and vice versa, see Fig. The reported results for the </?($) are 

obtained at the experimental values u/2n = Ur/27T = 7.78 GHz, Ip = 74 nA and 
the following parameters: decoherence rate F^/27r = 33 MHz, average qubit excitation 
gaps A/271 = 5.9 GHz with normal distribution and dispersion cta = 3.6%, effective 
qubit-cavity coupling g/27i = 1.1 MHz with small dispersion ag = 1% and N = 7 qubits 
in the ensemble. These parameters give best fit in the experimental regime at small N 
discussed in the next section 

Figure shows plots of the phase shift dependence on magnetic flux for N = 
20, 50,100, 250, 500 qubits in the system and fixed spread cta = 3.6% of the normal 
distribution of Aj. One can see that for N < 100 the phase shift fluctuates and depends 
on a particular realization of the qubit frequencies, while for N > 100 it appears rather 
smooth and regular, meaning that the array of qubits can be treated as a continuous 
system at these parameters. In Fig. we show results for different values of disorder 
with the following dispersions cta = 5, 7.5,10,15, 25%, keeping qubit number N = 250 
fixed. One can see from the Fig. |^that, with increasing disorder, the phase shift pictures 
start to reveal fluctuations. For the realistic spread ua = 10% the system remains far 
from the ensemble limit, even for this relatively large number of qubits. 

3.2. Experimental Regime, Small N 

Next, we compare experimental [1] (Fig. left panel) and numerical (Fig. right 
panel) results for the phase shift (p(<F). The resonator frequency a;/27r = ojr/27i = 7.78 
GHz and persistent current Ip = 74 nA are fixed. In our numerical method, we assume 
a normal distribution of the qubit excitation gaps and effective coupling. We selected a 
parameter distribution with dispersions cta = 3.6% and ag = 1% for a system containing 
N = 7 qubits with an average excitation gap A/27r = 5.9 GHz, which closely resembles 
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Figure 2: Numerical results for the phase shift (/?($) at qubit energy gap disorder 
cta = 3.6% of = 20, 50,100, 250, 500 qubits. 





Figure 3: Numerical results for the phase shift (p(‘F) of iV = 250 qubits using five 
different values of disorder cta = 5, 7.5,10,15, 25%. Parameters of the system: probe 
frequency a;/27r = Ur/^n = 7.78 GHz and relaxation rate k/2ti = 715 kHz, average 
excitation gap A/2tt = 5.9 GHz of qubits with normal distribution and dispersion 
(Ta = 3.6%, average coupling constant 5'/27r = 1.1 MHz with ag = 1%, decoherence 
F^/27r = 33 MHz, persistent current Ip = 74 uA. 

the experimental data. Subsequently, we htted the decoherence rate and effective qubit 
coupling and found F^/27r = 33 MHz and g/27r = 1.3 MHz, respectively. While the 
assumed spread is less than the experimentally expected one (cTa^ < 20 %) all found 
parameters are in correspondence with those quoted in jl], where A®^P/27r = 5.6 GHz, 
gcxp/ 2 p: — MHz, F®^P/27r = 55 MHz and = 8. The values for qubit number and 
dephasing reported in Ref. pQ were found under the assumption of identical qubits. Here 






















we show, that the experimental data can be reproduced fairly well under the assumption 
of randomly distributed qubit parameters. In the experiment it was observed, that 
an ensemble of = 8 qubits, resonantly interacting with the cavity mode and 

monitored over long time, spontaneously dissolved into two sub-ensembles of 4 qubits 
each, resulting into two jumps in Fig. (left panel). In our exact diagonalization 
procedure we do not find the formation of sub-ensembles but arrive at correspondence 
with the experimental curves if parameters of the metamaterial and disorder are close 
to the experimental values. 




<f> - (&0/2 (m< 3 >o) 


Figure 4: Experimental (left panel) and numerical (right panel) results for the photon 
phase shift </9($) as a function of external flux bias <h, calculated for the system of iV = 7 
qubits. 


4. Conclusions 

We theoretically studied the model of a flux qubit array coupled to a cavity, with 
disorder in qubit excitation frequencies. The system under consideration contains a 
finite number of qubits and operates in the intermediate regime where disorder range 
and decoherence rate exceed the effective qubit-cavity coupling. We calculated the 
photon Green function through the exact diagonalization of the Hamiltonian in the single 
excitation regime, assuming low power of external microwave driving. We found that the 
resonant phase shift of a transmitted probe signal shows quantitative correspondence 
with the experimental data [1]. Variations in phase shift characteristics at different 
values of the disorder and number of qubits in the system were presented. 
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